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Abstract 


We have developed a variational perturbation theory based on the Liouville- 
Neumann equation, which enables one to systematically compute the pertur¬ 
bative correction terms to the variationally determined wave functions of the 
time-dependent systems. We then apply the method to the time-independent 
anharmonic oscillator, and show that the results agree with those of other 
variational perturbation theories. We also show that the system has an inter¬ 
esting algebraic structure at the first order correction level. 
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I. INTRODUCTION 


Recently considerable efforts have been devoted to construct systematic procedures to 
improve the variational approximation method of quantum mechanics and quantum held 
theory [^]-|^. The so-called variational perturbation theory |l|J^ is based on the observation 
that the variationally determined approximate ground state of an interacting system is the 
exact ground state of a new system described by a quadratic (or, a soluble) Hamiltonian, 
and one uses the exact eigenstates of this new Hamiltonian to formulate a perturbation 
theory. Although this method is conceptually simple, the computation of the perturbation 
series is quite complicated for complex and especially time-dependent systems. The method 
of reference 0, on the other hand, is based on the background held method, and provides 
a fairly simple procedure to compute the perturbative correction terms to the variationally 
determined Gaussian ehective action. For time-independent systems, the latter method gives 
the same ehective potentials as those of the variational perturbation theory of ref. |]l|,|^. If 
one is interested in the explicit time evolution of a system of interest, however, one must 
construct explicit time-dependent wave function(al)s of the system, and the methods of ref. 
H are not of much help for this purpose. 

It is the purpose of this paper to develop another method of perturbative expansion 
around the variational approximation, which provides a systematic procedure to evaluate 
the correction terms to the variationally determined wave function(al)s of general time- 
dependent systems. 

To do this we use the technique of Lewis and Riesenfeld Q], which is originally devel¬ 
oped for solving the Schrodinger equation for a time-dependent harmonic oscillator. Let us 
consider the Schrodinger equation (unit h = 1), 

= ( 1 ) 

where the Hamiltonian H may have explicit time-dependence. The method of ref. ^ consists 
of introducing a non-degenerate Hermitian invariant operator, /, which satishes the Liouville- 
Neumann equation, 

iL/(()+ [/((), HI = 0. (2) 

It is the property of this invariant that the eigenvalues are time-independent, ^ = 0 [§]: 

(^) 

One can then show that the solution of the Schrodinger equation (Q) is given by 

|^(t)) = (4) 

n 

Therefore the task of hnding the solution of the Schrodinger equation is equivalent to hnding 
the invariant operators. 

In oder to develop a new perturbative method based on the Lewis-Riesenfeld technique, 
we will use a quantum mechanical system described by the Hamiltonian (|^ below as a 
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simple example. For this purpose it is useful to express the invariant operator, /, in terms 
of lowering and raising type operators a and ah 

m = E (6) 

n,r,s 

where a and also satisfy the Liouville-Neumann equation (2) for a new quadratic Hamil¬ 
tonian Hg and the commutation relation, 

[a,a^] = l, ( 6 ) 

and A is an appropriate expansion parameter. The quadratic Hamiltonian Hq will be de¬ 
termined such that it gives the variational approximation of the interacting system. Now 
finding the invariant operator /, which enables one to obtain the solution to the Schrodinger 
equation (|l|), amounts to finding the annihilation and creation operators which satisfy both 
the Liouville-Neumann equation (|^) for Hq and the normalization condition (^), and to 
finding the coefficients of Eq. (^. 

To obtain the variational approximation method for an interacting system we write a 
and af in terms of the dynamical variables of the system with an arbitrary complex function 
introduced, and express the Hamiltonian in terms of these a and ah We then truncate the 
Hamiltonian up to quadratic terms in a and for an approximation. The requirement that 
a and satisfy the Liouville-Neumann equation with the truncated quadratic Hamiltonian 
determines the equation for the variational approximation of the original system. We then 
treat the remaining part of the Hamiltonian perturbatively to determine the annihilation 
and creation operators of the interacting system as a perturbation series, which leads to the 
Liouville-Neumann approach to the variational perturbation theory. 

In the next section we establish the variational approximation method based on the 
Liouville-Neumann equation, by constructing a quadratic Hamiltonian which determines the 
variational approximation of the interacting system. In section HI we develop a systematic 
procedure to evaluate perturbative correction terms to the variationally determined wave 
function of the interacting system. In section IV we apply the method to the case of the 
eigenvalue problem for the time-independent anharmonic oscillator, and show that the result 
is equivalent to that obtained by other variational perturbation theories. We conclude with 
some discussions on our method in the last section. 


II. LIOUVILLE-NEUMANN APPROACH TO THE VARIATIONAL 

APPROXIMATION 


To develop the variational perturbative method based on the Liouville-Neumann equa¬ 
tion, we consider the time-dependent anharmonic oscillator described by the Hamiltonian, 


H 





( 7 ) 


where the mass is scaled to unity for simplicity, and ou and A are time dependent in general. 
It is convenient to express the time-dependent coupling constant A as 
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X{t) = \e{t) 


( 8 ) 


where A is the time-independent expansion parameter and e{t) carries all the time- 
dependence of the coupling constant A(t). For the time-independent systems we will set 
e{t) = 1. As explained in the previous section, the variationally approximated states of this 
system is the exact states of a quadratic variational Hamiltonian. To hnd such a quadratic 
Hamiltonian for the anharmonic oscillator (]^), we introduce a set of annihilation and creation 
operators 0, 

a = i(u*{t)p — u*{t)q'j , 

aX =—i(^u{t)p — u{t)q'j, (9) 

where u{t) is an arbitrary complex function satisfying the normalization condition, 

{u*u — iiu*) = i, (10) 

due to the commutation relation (|^). 

We now express the Hamiltonian (|^) in terms of a and of (|^), which, after normal 
ordering, becomes. 


H = Hg + XH', 


( 11 ) 


where 


and 


Hg = 


• Jf: • I 

U U + 


uj‘^u*u + ?)Xe{u*uY {a^a + 


- Xeiu*uY 

4 


2L 

1 r 
2- 


-1- + ?,Xe{u*u)u 


*2 


42 


-|- 3Xe{u*u)u^ 


( 12 ) 


h'-iJ: 


k=0 


.k 




(13) 


To hnd the quadratic Hamiltonian that determines the variational approximation, we require 
the annihilation and creation operators (^) to satisfy the Liouville-Neumann equation. 


d 

[a, ^g] = 0, 

+ [«^^G] = 0, 


(14) 


with the truncated Hamiltonian Hg of Eq. (|^. Eq. (|I^) requires the complex function 
u{t) of (^ and ([T0|) to satisfy the equation of motion. 


u + io'^u -|- 3Xe{u*u)u = 0. 


(15) 


4 










Note that Eq. ([T5|) is exactly the gap equation appearing in the Gaussian approximation 
where the trial wave functional is written as |P], 


^(g) = A^exp 


. U 2 

^- q 

L 2u* 


(16) 


To obtain solutions of the Schrodinger equation with Hq one needs to hnd the eigenstates 
of the invariant operator I = a^a, which consist of the state annihilated by the operator a, 


6|0)[o] ^ 0, 


(17) 


and the number states 




(18) 


In the coordinate representation the state, |0)[o], is represented by a Gaussian wave function, 

1 


'^[o],o{q,t) = ( 


27ru*{t)u{t) 


1/4 2 

exp t- — —q 


L 2u*{t) 


(19) 


Having found the eigenstates of the invariant operator, one can easily construct the general 
solution of the Schrodinger equation with hfc from Eq. (^). Note that Eq. (|T^ is the same as 
(0), which shows that our procedure indeed leads to the correct variational approximation 
of the system. The expectation value of the Hamiltonian with respect to this state is found 
to be 


3 

](0|.H|0)[o] = -(iL*u + uj‘^u*u + — -Xe{u*uy 


( 20 ) 


This method is also useful for solving the eigenvalue problem for the Hamiltonian of 
the time-independent anharmonic oscillator, in which case the parameters uj and A are now 
time-independent(e(f) = 1). In this case Eq. ([T3|) has a simple solution. 


u{t) = 


\/2^1q 


where Gg satishes the gap equation. 


+ 




3A 


( 21 ) 


2U 


G 


The state annihilated by the annihilation operator, a, is then given by, 

- /fUlVVinos" 




and the expectation value of the Hamiltonian, (pOD, reduces to 

Gg 3A 


[o](0|H'|0)[o] = 


16G 


2 ’ 


G 


( 22 ) 


(23) 


(24) 
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which coincides with the result of the Gaussian approximation of the anharmonic oscillator 
This indeed shows that our method reproduces the well-known results of the variational 
approximation. 

A few comments are in order. For the solution (|^) the terms proportional to and 
in Hg vanish due to Eq. (^2]), which shows that the variational Gaussian approxima¬ 


tion is the approximation of the interacting system by a simple harmonic oscillator with 
a new frequency, Qg- In spite of the apparent time-dependence of the complex function, 
u{t), the quantum ground state, is time-independent due to the choice of the solution 
(0), which minimizes the expectation value of the Hamiltonian, Eq. PD|). There also exist 
time-dependent nontrivial solutions to the equation (pAf ) for the time-independent anhar¬ 
monic oscillator, which correspond to the squeezed states in the case of the simple harmonic 
oscillator . Such solutions turn out to be convenient in studying the time evolution of the 
Universe at the early stage of the inflationary period |^. 

Finally we remark that in the case of the simple harmonic oscillator (A = 0), the solution 
(PI) gives the standard vacuum, and for a time-dependent harmonic oscillator, which can 
be regarded as the zero-mode of a minimal scalar held in a FRW background, it coincides 
with the Bunch-Davies vacuum [R]. 


III. PERTURBATIVE CORRECTION TO THE VARIATIONAL GAUSSIAN 
APPROXIMATION: SOLUTIONS OF THE SCHRODINGER EQUATION 

Having established the variational approximation method based on the Liouville- 
Neumann equation, we now proceed to develop a systematic method of computing the 
correction terms to the variationally determined solution of the Schrodinger equation. To 
do this we expand the annihilation and creation type operators, A and A\ which satisfy 
the Liouville-Neumann equation (^) with the total Hamiltonian, as a power series in the 
expansion parameter A: 


i = a + ^ ^ A'^H(„), 


n=l 

oo 


n=0 

oo 


it=at+ 


n=l 


n=0 


(25) 


where .B(o) = a, = a\ a and are the annihilation and creation type operators 
determined in the last section, and B^n) and are expanded as, 

G) = E( 26 ) 


with being the c-number function of time [Note that the operators B(^n) and Bj^^^ are 

dependent on e{t) through Substituting Eq. (P^f ) into the Liouville-Neumann equation 

(0), and by comparing the coefficients of A", one obtains the equation. 



in) 




(27) 
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where = 0. 

We note that any product of the operators that satisfy Liouville-Neumann equation also 
satishes the Liouville-Neumann equation. Thus, the operator satishes the equation, 

^(at’-a^) =^[at^a^iLG], (28) 


and Eq. (^Tl) reduces to 



(29) 


Note that the time derivative operator in the left hand side of Eq.(^) acts only on the 
coefficients This shows that we can End B(^n) by treating a and as if they are 

time-independent: one can integrate the equation (^^ to get 



■■■ (30) 

where the integrals are ordered from the inner bracket to the outer bracket, i.e., for example. 


B^2) = tdf f dt"[[a,H'{t")lH'{t')]. 

Jo Jo 

(31) 

The state annihilated by A, 


i|0) = 0, 

(32) 

can also be expanded as 


OO 

k=0 

(33) 

where 0)[o] is the variationally determined state, Eqs. (|^) and (|^, (5(o) 
an operator to be determined. Substituting (p3D into (|3^) we obtain 

= 1 and (5(fc) is 

OO 

i|0)= ^ A^+%Q(,)|0)[o] = 0. 

(34) 


/ c ,/=0 


For Eq. (|5^ to be satisfied, each coefficient of A” must vanish: 

k-\-l=n 

T. %Ow|0)|01=0, (35) 

k^l=Q 


which can be rewritten as 


7 


(36) 


n—1 

-®(0)Q(n) |0)[0] ®0(n)|0)[0] ^ ^ -^(n— fc)0(A:) |0) [0] • 

k=0 

By operating on both sides of Eq. (^), one finally obtains the relation, 

n—1 

Q(n)|0)[0] = ^ -S(n-fe)Q(fc)|0)[0]! (37) 

k=0 


where N is the nnmber operator, did. Note that the sum on k in the right hand side of Eq. 
(0) runs only up to (n — 1), and thus we can construct the operator (5(n) by applying Eq. 


(^) iteratively, starting from the hrst non-trivial contribution. Substituting (^) into (^) 
we obtain the state annihilated by A as a perturbation series, from which solutions of the 
Schrodinger equation can be constructed from Eq. ( 0 ). 


IV. EIGENVALUE PROBLEM OF THE TIME-INDEPENDENT ANHARMONIC 

OSCILLATOR 

To illustrate the effectiveness of our approximation method we now compute the ground 
state of the time-independent anharmonic oscillator described by the Hamiltonian (|^) with 
time-independent parameters uj and A (e(f) = 1). 

For the time-independent systems we are interested in the ground state of the system, 
which is an eigenstate of the Hamiltonian of the system. In the introduction we have shown 
that the solution of the Schrodinger equation (]^) can be written as a linear combination of 
the eigenstates, |e„,f), of the invariant operator J, which do not necessarily coincide with 
the eigenstates of the Hamiltonian. If the eigenstates are time-independent, ■^\en,t) = 0, 
however, one can easily show from Eqs. and (|^) that such states are the simultaneous 
eigenstates of I and H: 

[Hj]\en,t) = 0, (38) 

and 

( 39 ) 

from the non-degeneracy of I. It follows from this that the state annihilated by A, Eq. (P3[), 
is indeed the ground state of the time-independent system if it is time-independent. 

We have shown in section H that the state annihilated by a for the time-independent 
anharmonic oscillator, Eq. (p3D, is time-independent. This fact guarantees that the state 
(0) is indeed the variationally approximated ground state and the expectation value of H 
with respect to this state is the approximate ground state energy of the system. And the 
state annihilated by A, Eq. (^), for the time-independent system is the perturbatively 
corrected ground state if it is time-independent. 

The ground state energy of the time-independent anharmonic oscillator can be obtained 
from the ground state (0), 

OO 

E, = i: ( 40 ) 

n=0 
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where 


_ 3A 

Eg,m - ^ 

Eg,{i) = [o](0|if'|0)[o] = 0 (41) 

n—1 

Eg,{n) [0] (0|-^ Q(n—1) |0)[O] -^g,(fc) [0] (01Q(n—fc) 10) [0]; fl ^ 2. 

k=0 


Note that, due to the property of the variational perturbation theory, the hrst order correc¬ 
tion to the energy eigenvalue vanishes. 

To show the time-independence of the quantum ground state, and to compare with the 
standard method, it is convenient to write the operators in Eq. as 


a = a'*' = e 


(42) 


where ho and hg are the time-indepen dent, standard annihilation and creation operators. It 
is manifest that the operators in Eq. (|42|) satisfy the Liouville-Neumann equation (^ with 
the truncated Hamiltonian which can be written as 


Hg = - ) = 


11^(00 «o + 


The remaining perturbation Hamiltonian becomes 


H' = 


{ 1 ^ 0 ^ + ho'^ho -I- -ho^ho + agho + 




(43) 


(44) 


From now on we shall work with both the time-dependent operators a, d) and time- 
independent ones hg, hg, whenever it is more convenient. 

We now want to compute the operators A and A^ and the vacuum state. It is more 
convenient to work with the time-independent, standard operators, in terms of which Eq. 
(^) reads 


E 


^-i{r-s)nct 


E 

di 



^ fr ^ s 
% % 


^E 



(45) 


It is straightforward to compute the hrst-order correction term of the annihilation operator 




hg -|- 


A 


(2Hg)H4Hg 


-t3 , 

Og -|- 


hg hg -|- (3ft -|- Cq^Qi^CL, 


2n 


G 


2U 


G 


(46) 


where Cg is an integration constant. is the Hermitian conjugate of H[i]. We thus hnd, 
from Eq. (^3]), the hrst order vacuum state. 


|0)[i] — |0)[o] + AQ(i)|0)[o] 


( 47 ) 
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and the corresponding vacnnm energy, 




3A 


3A2 


leu 


G 




+ 0{X^ 


(48) 


which is the same as that obtained by other variational pertnrbation theories Note 

that the hrst order vacnnm state |0)[i] is time-independent as it shonld be. 

We now consider an interesting algebraic strnctnre that the hrst order variational per¬ 
tnrbation of the a n h armonic oscillator exhibits. With the specihc choice of the integration 
constant, cq = |, we hnd that to the order of A the fnll Hamiltonian can be expressed as. 


O A 

H = + ^ - + 0(A2), (49) 

and the commutation relation between A and A'^ becomes, 

SA 

[H[i], HJ'ij] = 1 -h + O(A^) 

G 
Q A 

= 1 + C>(A^). (50) 

The algebraic strnctnre of Eq. (|49|) for the Hamiltonian explains why the state |0)[i] annihi¬ 
lated by H[i] is the vacnnm state for the fnll Hamiltonian np to the order of A. It shonld be 
noted that when the invariant operators and the fnll Hamiltonian are truncated to the order 
of A^, the truncated Hamiltonian (|4^) has a factorized form and reminds us of a possible 
underlying supersymmetry and a deformed commutator (^). The deformation parameter is 
greater than one, so the deformed algebra does not interpolate between bosons and fermions. 
However, the deformed algebra can be used conveniently to hnd the excited quantum states. 


V. DISCUSSIONS AND CONCLUSION 

We have developed a variational perturbation theory based on the Liouville-Neumann 
equation, which enables one to compute the order-by-order correction terms to the varia- 
tionally approximated wave functions of the time-dependent systems. The method is based 
on the observation that the variationally determined approximate eigenstates of the Hamil¬ 
tonian of an interacting system are the exact eigenstates of a new quadratic (or, a soluble) 
Hamiltonian, and on the fact that the solution of the Schrodinger equation for a time- 
dependent harmonic oscillator can be represented by eigenstates of an invariant operator 
satisfying the Liouville-Neumann equation (^. We thus write the Hamiltonian of an inter¬ 
acting system as a function of annihilation and creation type operators a and a+, which are 
dehned as linear functions of dynamical variables with an arbitrary complex function u{t) in¬ 
troduced. We then truncate the Hamiltonian up to quadratic terms in a and a^, and require 
a and a"*" to satisfy the Liouville-Neumann equation with this truncated Hamiltonian Hq- 
This requires the function u{t) to satisfy a gap equation, which is exactly the same as that 
appearing in the Gaussian approximation method as shown in section H. This constitutes 
the Liouville-Neumann approach to the variational Gaussian approximation method. 
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To develop a perturbation theory based on this variational approximation, we expand the 
annihilation and creation type operators A and A'^, for the full interacting Hamiltonian, as 
power series in the coupling constant A. By requiring the operators A and A"*" to satisfy the 
Liouville-Neumann equation with the total Hamiltonian, these operators A and A'^ can be 
determined to each order in the coupling constant A. By constructing the state annihilated 
by the operator A, we obtain the quantum state satisfying the Schrodinger equations to each 
order in A. Since this method enables one to construct wave functions of time-dependent 
systems, it can be conveniently used to study the time-evolution of time-dependent systems 
such as time-evolution of scalar held in the expanding Universe. 

To hnd the eigenstates of time-independent Hamiltonian, one needs one more requirement 
to the procedure explained above. For the eigenstates of the invariant operator I to be the 
simultaneous eigenstates of the Hamiltonian, the eigenstates must be time-independent as 
shown in the last section. We have used this method to compute the ground state of the 
time-independent anharmonic oscillator to the hrst order in A, and have shown that the 
result agrees with that of other variational perturbation theories. To this order of the 
perturbation the Hamiltonian is found to be factorized as in the simple harmonic oscillator 
with the annihilation and creation operator A and A'^, which satisfy the deformed algebra 
(0)- This interesting algebraic structure can be used to construct all the excited states of 
the system. This possibility will be explored elsewhere [^. 

Although we have considered the quantum mechanical anharmonic oscillator for simplic¬ 
ity, this method can easily be generalized to the case of quantum held theory. 
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